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I. INTRODUCTION 

Over the last decade or so there has been a great revival of interest in models featuring extra dimensions of space. Kaluza- 
Klein theory, originally proposed early in the 20th century, is the forerunner of modern universal extra-dimension models [HIH]. 
These theories hypothesise additional small extra dimensions of compact topology, with all degrees of freedom unrestricted in 
their propagation. An alternative scenario sees some fields localised to a submanifold or brane, while other fields propagate 
throughout the bulk. In particular, Arkani-Hamed, Dimopoulos and Dvali pointed out that compact extra dimensions may be 
rather large in characteristic scale if only gravity is able to propagate in the bulk 

At about the same time, two seminal papers by Randall and Sundrum, referred to hereinafter as RSI and RS2, investigated 
brane-world models with non-trivial spacetime curvature in the bulk (warped metrics). RS 1 is a model incorporating one compact 
extra dimension with a topological identification [y] e /Z2 and a warped-metric geometry [4]. Two branes, denoted Planck 
and TeV, are placed at the orbifold fixed points [y] = and [y] = FcTt, respectively. The original setup had all standard-model 
fields confined to the TeV brane, while the graviton was dynamically-localised to the Planck brane. The metric warp factor 
aggregates normalisation conditions so that the physical mass scale, on the TeV brane, is related by an exponential hierarchy to 
the free parameters in the Lagrangian. Thus it is possible to have the Lagrangian parameters at a high scale while maintaining 
comparatively small standard model masses, thereby solving the gauge hierarchy problem. RS2 is obtained, effectively, by taking 
the TeV brane infinitely far from the Planck brane and reflecting the resulting spacetime about the Planck brane |5]. The result 
is a model with an infinite extra dimension but still with effective four-dimensional gravity dynamically induced on the brane. 
The standard model fields must now be taken to live on the Planck brane. [See references L6J,L7J,L8J for other foundational work 
on extensions to the standard model using extra dimensions.] 

RS2 is of considerable interest because the extra dimension is treated analogously to the usual three dimensions of space, 
unlike the compact extra dimension paradigm. However, pure RS2 leaves open the question of what the brane is (D-brane of 
string theory?) and why non-gravitational fields are localised to it. Also, the extra dimension is not exactly on the same footing 
as the other dimensions, because the placement of the fundamental brane explicitly breaks translational invariance along the 
extra dimension. 

Independently of the above, it had been proposed that our universe might be a domain wall or other topological soliton [9] . See 
also ifloll . lfTTll . In the codimension 1 case of a domain-wall, the setup is quite similar to RS2. Such a universe is brane-like, but 
it has finite thickness and the origin of the brane is specified. Brane formation is now an instance of the spontaneous rather than 
the explicit breakdown of translational invariance. Finally, one can envisage that all fields are localised for dynamical reasons to 
the brane, not just the graviton as in pure RS2. What, exactly, would such a theory look like? 

Recently, one of the authors (RRV), proposed, together with Davies and George, a specific theory (to be called the DGV 
model hereinafter) that may realise the dynamical localisation of an SU(3)c x SU(2)iv x U(l)y gauge theory plus gravity to a 
solitonic domain-wall brane with one warped extra dimension 1 12]. The DGV paper speculates that Dvali-Shifman gauge field 
localisation (to be reviewed in a later section) is implemented by confining SU(5) bulk gauge dynamics. The purpose of this 
paper is to extend this model by enlarging the bulk gauge group to SO(IO). This allows the domain-wall localised gauge theory 
to be either SU(5) x U(l)x or the left-right symmetric model, depending on the choice of Higgs potential. 

In the next section, we discuss the qualitative features of this kind of domain-wall brane model, as a warm up for the detailed 
constructions presented in sections [TVl and [Vl Section |III] critically reviews the Dvali-Shifman mechanism, while our final two 
sections contain a discussion of salient points for our model and conclusion. 

II. THE DVALI-SHIFMAN DOMAIN- WALL BRANE SETUP 

The most basic structure in domain-wall brane models is a kink configuration for a scalar field. The prototype kink is supplied 
by a Z2 invariant (p'^ potential for a real scalar field 0. The Euler-Lagrange equations possess a solitonic kink solution, 

4>{xf^,y) - vtanh(m3;), (1) 

where v is the Z2-breaking vacuum expectation value, while m is the inverse width of the kink. Both parameters are given 
through the Higgs potential. The extra dimension is described by the coordinate y. The energy density profile for this solution 
is localised around the zero of the kink (y - 0). We say the scalar field has condensed to form a domain wall. The region 
away from the 3 H-1 -dimensional domain wall is called the bulk. The configuration is topologically stable because it asymptotes 
to vacuum expectation values that spontaneously break a discrete symmetry, and are thus disconnected from each other in the 
vacuum manifold. 

The prototype domain wall must be embedded in a considerably richer theoretical structure in order to produce a model that 
might have realistic phenomenology. In particular, the requirement that massless gauge bosons be dynamically localised to the 
wall is most plausibly met by invoking the Dvali-Shifman mechanism f\3\. This proposal will be critically reviewed in the next 
section. It requires the unbroken gauge group, H, inside the wall to be a subgroup of that in the bulk G. The gauge bosons of H 
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are then thought to be localised, provided that the bulk gauge theory is in a confinement phase. 

The DGV model uses G - SU(5) and H - SU(3)c x SU(2)vt' x U(l)y, thus potentially realising a domain-wall localised 
standard model. It is also a novel reinterpretation of SU(5) grand unification. As with standard SU(5) theories, the extension to 
SO(IO) immediately suggests itself because all the standard model fermions of a given family may be assembled into a single 
irreducible representation, the 16. In this paper, we thus start with G - SO(IO). We also need a discrete symmetry outside of 
SO(IO) for topological-stability reasons; we choose Z2 for simplicity and economy.' 

We therefore consider an SO(IO) x Z2 gauge-invariant Lagrangian, eventually to be coupled to a Randall-Sundrum warped 
4-1- 1 -dimensional metric in order to dynamically localise the graviton. Group-theoretically extrapolating the original argument 
presented by Dvali and Shifman, an SO(IO) singlet scalar field in conjunction with an SO(IO) adjoint Higgs is used to dynami- 
cally generate the required domain wall. The singlet Higgs takes on a kink configuration, while the adjoint Higgs configuration 
is nonzero inside the wall, spontaneously breaking SO(IO), and asympototes to zero at \y\= 00 thus restoring SO(IO) in the bulk. 
For quartic potentials, we shall show that the unbroken group inside the wall is SU(5) x U(l)x, while the extension to sixth order 
permits the further breaking to the left-right group SU(3)c x SU(2)i x SU(2)r x U(1)b_£. If the Dvali-Shifman mechanism is 
operative, then the gauge bosons of these subgroups are dynamically localised. _ 

Chiral 3-1-1 -dimensional fermions are localised by their interactions with the background Higgs field configurations [9]. Each 
fermion is confined around a 3-1-1 -dimensional hypersurface parameterised by certain values for the bulk coordinate y. These 
values differ according to fermion species, because the adjoint Higgs configuration "splits" their localisation points. 

For completeness we discuss subsequently breaking SU(5) x U(l)x and, respectively, SU(3)c x SU(2)l x SU(2)j; x \J(1)b-l 
down to the standard model on the domain wall brane. This will require additional Higgs fields which must condense inside the 
domain wall and some fine tuning conditions to induce the required gauge hierarchies. 

An important advantage of this type of 4+ 1 -dimensional grand-unified theory is that the usual tree level mass relations between 
the quark and lepton masses no longer appear This is precisely because the fermions are split along the extra dimension. Each 
fermion's 3+1 -dimensional mass scale depends on overlap integrals of the fermion's bulk profile with the extra-dimensional 
profile functions of the additional Higgs fields which we introduced to break the symmetry down to the standard model gauge 
group. For the SU(5) x U(l)x case, the bulk profiles for the fermions depend on their distinct U(l)x charges, thus these overlap 
integrals contribute factors to the 3-1-1 -dimensional mass parameters which are different for fermions in different SU(5) x U(l)x 
representations. Hence the tree level mass relations are reduced to the SU(5)xU(l)x subset of the normal SO(IO) mass relations. 
A similar effect occurs for the left-right symmetric alternative model. 

Before discussing the technical details of the above models, we need to pause and talk briefly about using the Dvali-Shifman 
mechanism to confine gauge fields to a codimension 1 domain- wall brane in 4-1-1 -dimensions. We shall use this opportunity to 
offer a brief heuristic review of the mechanism. 



m. DVALI-SHIFMAN GAUGE FIELD LOCALISATION 

It is straightforward to dynamically localise both chiral fermions and other scalar fields to a domain-wall background. But 
the localisation of massless gauge fields is a much greater challenge. The simple analogue of the fermion and scalar localisation 
mechanisms can produce localised massive spin-1 fields, but this is not what we require, both because photons and gluons 
are massless, and because gauge coupling constant universality is lost due to differing overlap integrals involving the extra- 
dimensional profile functions. 

A very different mechanism is required, and the most well-known proposal is conceptually diflferent from the purely classical 
physics one may use to trap fermions and scalars. DvaU and Shifman consider a toy model based on an SU(2) Yang-Mills theory 
in 3+1-dimensions breaking to U(l) on a 2-1-1 -dimensional domain wall. The bulk is presumed to be in confinement phase with 
a strong coupling regime in the infrared. This creates a mass gap of order the bulk confinement scale Aconf[su(2)] in the vector 
boson spectrum, so away from the wall the gauge fields form bound states with mass > Aconf [SU(2)] ■ 

On the domain wall the photon is free, however in order to propagate in the direction transverse to the wall it must become 
conglomerated into a massive SU(2) glueball. The energy cost associated with moving off the wall is untenable for the massless 
U(l) gauge boson which effectively becomes trapped on the 2-1-1 -dimensional topological defect. 

In the 't Hooft-Mandelstam flTlflS^ dual superconductivity picture |19], electric field lines from test charges on the brane 
spread out like normal Coulomb fields along the brane directions. However, when they meet the bulk they get repelled by the 
dual analogue of the Meissner effect. Thus any two charges located on the brane will interact via a 2-1-1 -dimensional Coulomb 
force at distances much greater than the thickness of the wall. 



An earlier attempt at SO(IO) domain-wall brane models can be found in fT4ll . based on the clash of symmetries idea fT^ . However, these schemes are unable 
to produce phenomenologically acceptable fermion localisation fl^ . 



4 



If the test charge is located in the bulk then the field lines are no longer at liberty to spread out. Instead they form a flux 
string which tunnels through to the brane and expels the field onto the domain wall. This is the dual superconductor effect 
corresponding to Abrikosov vortex formation. Hence, regardless of the bulk localisation profiles of the fermions, the effective 
field on the domain wall will exhibit the same 2+1 -dimensional Coulomb distribution. 

The Dvali-Shifman argument can easily be extended to deal with a more general gauge group G breaking to some subgroup 
H on the domain wall. Under these circumstances it is assumed that the confinement scale on the brane, Aconf[//], is significantly 
lower than the bulk scale, Aconf[G] • 

However, extending the model to cope with a larger number of dimensions is not easy. The problem is that non-Abelian gauge 
theories are not renormalisable in more than 3+1-dimensions. 

3+1 -dimensional non-perturbative lattice simulations provide qualified support for the hypothesis that the Dvali-Shifman 
mechanism can trap gauge fields on a codimension-1 brane |20]. But to our knowledge no one has tackled a 4+1 -dimensional 
simulation. At this stage the Dvali-Shifman mechanism remains an intriguing conjecture in 4+1-dimensions. 

A necessary condition for the Dvali-Shifman mechanism is confinement in the bulk. It is encouraging that lattice gauge 
theory simulations for SU(2) pure Yang-Mills theory in 4+1-dimensions demonstrate a first order phase transition at finite lattice 
spacing as a function of the gauge coupling constant: for coupling strengths above a critical value, the theory appears to be 
confining [21,1 . This result has been extended to pure SU(5) Yang Mills gauge theories in [22]. 

However no conclusion can be drawn about the continuum limit since our non-renormalisable theory has problems with point- 
like interactions at high energies. But this may not actually be a problem for our application. All field-theoretic brane models 
are necessarily low-energy effective theories due to their non-renormalisability, so they are implicitly defined with an ultraviolet 
cutoff. The finite lattice spacing in the simulations is also an ultraviolet cutoff, so it appears sensible to use the lattice gauge 
results to conclude that confinement in 4+1 -dimensional effective gauge theories can exist. 

While this argumentation provides encouragement to pursue models based on the Dvali-Shifman idea, it does not rigorously 
establish that any specific model works. For one thing, no such model is a pure Yang-Mills theory. Secondly, for each candidate 
theory one would need to compute via lattice simulations the critical gauge coupling constant as a function of the cutoff. Since 
the critical coupling constant in 4+1 dimensions is itself of nonzero mass dimension, one would then need to check that this 
value is compatible with other scales in the problem, such as the inverse width of the domain wall and the effective grand unified 
symmetry breaking scale inside the wall. In particular, we need the SO(IO) breaking scale to be higher than the bulk confinement 
scale. 

We find it interesting that the assumed validity of the Dvali-Shifman mechanism allows plausible candidate domain-wall 
localised gauge theories to be readily constructed. We believe this provides good motivation to pursue more in-depth studies of 
the mechanism via lattice gauge theory. 

IV. THE SO(10)-»SU(5)xU(l)x MODEL 

Our goal is an SU(5) x U(l)x gauge theory on the 3+1-brane with a full SO(IO) unified theory in the bulk. For simplicity we 
shall start with a Minkowski flat action. We then compile our theory using the following algorithm: 

• construct a domain wall using a scalar SO(IO) singlet; 

• use an adjoint Higgs to break SO(IO) to SU(5) x U(l)x on the domain wall and invoke Dvali-Shifman gauge field locali- 
sation; 

• confine zero mode chiral fermions; 

• add a Randall-Sundrum warped metric; 

• discuss breaking SU(5) x U(l)x to the standard model on the domain wall and electroweak symmetry breaking. 

It is helpful to establish a consistent notation before we start. The indices M and N run from 1 to 5, while fi, V denote 
coordinate indices belonging to the subspace 1 to 4. We distinguish the coordinate - y. Our convention for the metric 
signature is diag(+, -). Elements of the set {0-1,0-2,0-2} are used to label the 2x2 Pauli matrices, or 3 dimensional 
Clifford algebra. 
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A. Domain wall construction and gauge field localisation 

We consider a Higgs sector constituted by an SO(IO) Higgs singlet <p ~ 1 and an adjoint Higgs X ~ 45. We impose the 
discrete Z2 symmetry, y -y, — > - 0, A" — > -A". The Higgs potential is 

y_^_^ ^ - yTrA^ + ^(TrX^f + ^TrX' + '^cp^TvX^ + ^(0^ _ v^)^, (2) 

where we have truncated our potential at fourth order. We argue that the general expression for Vx,,/, can be expanded as a 
polynomial in the fields X and (/> where higher order terms are suppressed in an effective low energy theory by their dimensionful 
coupling constants. We include 4th order terms because orthogonality of the SO(IO) generators remonstrates that truncating 
our potential at second order would have introduced an accidental symmetry whereby each of the 45 SO(IO) adjoint-Higgs 
field components could be transformed independently under SO(IO). Since our theory is non-renormalizable in 4+1 -dimensions 
and is therefore considered to be effective only up to an ultraviolet cut off we do not expect pathologies from the dimensionful 
coupling constants. The adjoint Higgs multiplet is represented by the 10 x 10 anitsymmetric matrix X - 2/ TjXi where the T, 
are the generators of the fundamental representation normalized so the Tr (TiTj^ = j6ij. In a certain regime of parameter space 
the vacuum manifold is (X = Q,<p = +V). Solitonic solutions to the Euler-Lagrange equations must obey Dirichlet boundary 
conditions, where we require both Higgs fields to asymptotically approach vacuum configurations. We choose the boundary 
condition to be (0, -V) at y = -00 and (0, V) at y = +00. 

Since the vacuum manifold is not connected, the relevant homotopy group ;ro(SO(10) X Z2/SO(10)) is non-trivial. Our 
topological boundary conditions exploit this. Any solution to the Euler-Lagrange equations satisfying these boundary conditions 
cannot spontaneously evolve into a solution from another topological sector, irrespective of relative energy densities. Id est, 
any map interpolating between the two disconnected pieces cannot be deformed continuously to give the trivial map from S'*' to 
SO(10)xZ2/SO(10). 

Under these conditions, the stable solution to the Euler-Lagrange equations is the minimum energy density solution belonging 
to this topological sector. Like Dvali and Shifman, we identify stable solutions by checking the dynamical evolution of pertur- 
bative linear modes. We also checked our results numerically to account for higher order effects. The energy density of each 
solution is dependent on the coupling constants, hence different solutions are stable in different coupling constant regimes. 

We find that A2 = is a bifurcation point for the manifold of solutions to the Euler-Lagrange equations. When A2 > the 
SO(IO) adjoint field will condense about y = 0, so that the nonzero components of X arrange for SO(IO) — > SU(5) x U(l)x on 
the domain wall. When A2 < our analysis indicates that to first order the X field components will try to adopt a configuration 
which breaks SO(IO) — > U(l)^ on the domain wall. A more detailed discussion is given in the Appendix dVIIIb . 

Let us take a closer look at the A2 > scenario. Solutions to the Euler-Lagrange equations persist for a wide range of 
parameter values. Purely for the sake of convenience, however, we choose to focus on a concrete analytic solution that exists 
provided the parameters obey 

(4//^ - 2k +A\/^)(20Ai + lAi) - 5k{Ai? - kV^) = 0. (3) 

We emphasise that this is not a fine-tuning condition. Rather, it defines a special slice through parameter space that happens to 
admit an analytical solution. 

Under these circumstances the Euler-Lagrange equations have a stable solution of the form, 

(p(y) - V tanh(my), Xi - Asech(my), (4) 

with = "^To/h+liT ' ~ ~2;U^ + kV^ and Xi is associated with the U(1)a' generator which we embed inside SO(IO) according 
to Ti = -^diag(cr2, 0-2, cr2, cr2, 1x2), where 0-2 is the second Pauli matrix. Aside from having a closed form, this solution 
is convenient since the dynamical equations for the first order terms in a perturbative expansion about this solution can be 
transformed into hypergeometric differential equations. We are able to solve the equations exactly and show that all perturbative 
modes are oscillatory. So, the solution is stable against further condensation of the remaining 44 X field components. 

As promised, this solution spontaneously breaks the SO(IO) gauge symmetry down to SU(5) x U(l)x on the domain wall. So 
if we assume the Dvali-Shifman mechanism to be effective, that is assume the 4+1 -dimensional bulk SO(IO) gauge theory is in 
a confining phase which is valid up to some ultraviolet cutoff Auv, then SU(5) x U(l)x gauge field locahsation will follow as a 
consequence of the background Higgs field configuration. 

Having established the existence of a stable Dvali-Shifman domain wall solution, we now turn our attention to the localisation 
of fermions. 
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B. Localising Fermions 

We extend the Dirac algebra to encompass 4+1 -dimensions by identifying F*' = (y'', -ij^), where {r'^,r*') = 27;*"^. The 
fermions are contained in a *P ~ 16 representation of SO(IO), and have the Lagrangian 

^Yukawa = /W^5m*P " h;^^'T"T''Xab'i' " h^Cp^^ , (5) 

where Xji^ is the jjcth entry of the 10x10 antisymmetric matrix used to represent X and W is a member of the 10 dimensional 
Clifford algebra 12311 . Specifically the Clifford algebra of dimension n is constituted by any n matrices satisfying {W, t*) = 
26j,k V;;A:e{l,...n). 

We choose to adopt the explicit form 

T^''-^ = 1 (g) ■ ■ ■ ® 1 (g) CTi (g) 0-3 (g) ■ ■ ■ ® 0-3, (6) 
T^* = 1 ■ ■ ■®1 ® 0-2 ® 0-3 (g) ■ ■ ■ (g) 0-3, (7) 

k-l 5-k 

where it is understood that the index k runs over { 1 , . . . 5). 

To maintain the discrete reflection symmetry we now require — » /r^*F in conjunction with y -y, (p -(f> and X — > -X. 

Each spinor ^'i(x,y), where the index / denotes the different irreducible SU(5) x U(l)x components of the 16, can be decom- 
posed in terms of a complete set of simultaneous eigenfunctions i//„Li{x) and il/nmix) for the 3 + 1-dimensional Dirac Hamiltonian 
and (to remove the degeneracy) the 3+ 1-dimensional chirality operator 7^. The basis states which appear with non zero coeffi- 
cients in our expansion will be the 3+ 1-dimensional chiral zero mode, if/oi or ip\^R depending on the case, plus a finite number of 
discrete massive modes, as well as a continuum of massive modes, as per 

y) = ^{fnu(y)Mx) + fnRi(y)^nRi{x)), (8) 

where the sum is understood to include an integration over the continuum. This is a dimensional reduction or generalised Kaluza- 
Klein procedure whereby the dependence on the y coordinate of "^lix, y) is subsumed into a complete set of mode functions or 
profiles fn,L/R,i, and the 4-1- 1-dimensional field is redescribed as an infinite tower of 3-1- 1-dimensional fields. The discrete mode 
functions are square-integrable, while the modes from the continuum are delta-function normalisable. 

We want the 4h- 1-dimensional Dirac fermions in the 16 of SO(IO) to supply 16 left-handed 3 H- 1-dimensional zero modes to 
be identified with a family of quarks and leptons. To display these modes, it suffices to truncate the full mode decomposition 
to just the chiral zero mode term fQu(y)^Qu{x)- This ansatz is now substituted into 4+1-d Dirac equation with the background 
domain-wall configuration playing the role of the mass term: 

= iY^'du'Vix, y) - hxT'^T^Xatiymx, y) - h^4>{ymx, y). (9) 

The separation of variables allows us to isolate the dependence on the y-coordinate, after requiring ffiat each of the if/Qu be left 
chiral and obey the usual massless 3+ 1-dimensional Dirac equation. Fortunately, foxX - Asech(OT3')ri the matrix T^T^Xatiy) is 
diagonal. We label the diagonal entry belonging to row / and column / by {j^'T''Xab{y))m where there is no intended sum over 
the index /. The differential equations decouple and the bulk localisation profiles for the fermions are then easily found to be 

where A^, is a normalisation constant. We choose the signs of the coupling constants, and h^, so that the integrand cuts ffie axis 
with a positive slope. The asymptotic behaviour of the kink implements the localisation of each fermion to a 3-1- 1-dimensional 
hyperplane coplanar with the zero of its respective intergrand, which occurs at y = yo- Due to coupling to the SO(IO) adjoint 
Higgs, fermions belonging to different SU(5) representations are localised around different parallel hyperplanes in ffie bulk, see 
Figure [T] 

These 3 H- 1-dimensional localised massless fermions are our SU(5) brane world matter. In flat spacetime ffiere is a mass gap 
between the zero mode and ffie lowest allowed massive 3+ 1-dimensional mode. Hence, at low energies, interactions between 
3 -Hi -dimensional zero mass modes can only produce other zero mass fermions l'24'l. This, in conjunction with the Dvali- 
Shifman gauge boson localisation conjecture, produces a candidate low-energy theory with dimensional reduction down to 3+1 
dimensions. 

We need to explain the physical interpretation we give to individual terms in our mode decomposition. Equation dS) is simply 
a convenient way of projecting the 4+1-d field onto a collection of 3+1-d coefficients iffn.L\Rj(x) by resolving the y-dependence 
of the field in terms of a complete set of orthonormal functions f„^L\R,iiy) which span the 1-dimensional rigged Hilbert space. We 
have done this because it allows us to easily see how 3+ 1-dimensional states with definite mass feel the classical background. It 
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FIG. 1 : The graph displays the integrand in the exponent of equation l IlQt which gives the extra dimensional profiles of the left handed zero 
modes associated with first generation fermions in the SU(5) x U(l)x brane world model. The graph shown here corresponds to the parameter 
choices = = I and the flat space background solution to the Euler-Lagrange equations Xi = Asech(m>'), = Vtanh(m>') with A = V = 1. 
3+1-dimensional left chiral fermions will be confined to the hyperplane corresponding to the zero of their integrand. 



is necessary to treat the different 3+1-dimensional eigenfunctions of independently because we have associated with dy. 
However mathematically equation ^ is simply a convenient way of writing the field. 

Since the physics is basis independent, we can choose to expand a general field H(x,y) = I.„gn(y)^„(x) (x = x'^) in terms of 
any complete set of states [gniy)}- For example: if the 'B.{x,y) is an SO(IO) gauge singlet then without loss of generality, we 
may conveniently choose the {^,,(3')) so that the coefficients {^„(x)) are solutions to the massive 3+1-dimensional Klein-Gordon 
equation. In this case the profile functions g„(y) appearing in the above expansion correspond to physically states of the effective 
3+1-dimensional theory which propagate as free particles in the confining bulk SO(IO) gauge theory. In our case the t//n.L\RAx) 
are charged under U(l)x and will be forced, by a confining bulk, to propagate as constituent particles of a bound SO(IO) gauge 
singlet. This is conceptually akin to writing the QCD lagrangian in terms of quarks and gluons with the implicit understanding 
that the propagating states are hadrons. Since SO(IO) is broken inside the wall, the 3+1-d fields behave in a manner consistent 
with the standard model gauge interactions. 

We argue that because our zero mode profile functions are sharply peaked around the domain wall and the SO(IO) confinement 
dynamics are suppressed here, our classical localisation profiles give us a reasonable first order approximation for the landscape 
of these low energy particles. See fl2ll . ll22ll for more details. 



C. Adding Warped Gravity 

The last vital component of our model is warped gravity. We search for solutions to the Einstein-Klein-Gordon equations for 
the action 

5 = J d'x ^{-2M'R -h + T-Vx^+ Xvukawa), (11) 

where G is the determinant of the five-dimensional metric tensor, M is the 5-d fundamental gravitational mass scale, R is the 
Ricci scalar, A is the bulk cosmological constant and T is simply the gauge field and Higgs boson kinetic terms. We would like 
a solution with 3+1 -dimensional Minkowski space and zero mode gravitons localised on our field theoretical brane. To achieve a 
Minkowski metric on the brane we impose a fine tuning condition on the bulk cosmological constant. The necessity of balancing 
the bulk cosmological constant against the brane tension is a feature of Randall-Sundrum like models. At the same time we need 
the qualitative forms of both X and to be similar to the flat-space case. Substituting a Randall-Sundrum warped metric ansatz, 

ds^ = g-pOO/fiM' Tj^^dx^'dx" -dy^, (12) 

into the Einstein-Klein-Gordon equations produces four coupled second order dififerential equations which exhibit appropriate 
solutions for a wide range of parameter values. By choosing, as before purely for convenience, our parameters to lie on the 
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manifold, 

em^M^ + V^m^ - 3M^A +3M^kV^ = 0, 

-Sot^M^ - - 6MV + 3M^Vv2 = 0, (13) 

20Ai +2A2- 10k + 5A = 0, 

we find an analytic solution of the form: 

p(y) - ^'in(cosh[my]), (f)(y) - W tanh[my]. A": = 2 V sech[my]. (14) 

This solution is consistent with all the localisation properties mentioned above. It provides the same Higgs configuration which 
we postulate induces Dvali-Shifman SU(5) x U(l)x gauge field localisation on the brane. 

In changing from flat space to a warped metric we must include vielbeins in the fermion action so that fields written at each 
point in terms of a local Lorentz coordinate system in the tangent space transform correctly under general coordinate transfor- 
mations [25J. We have also included the spin connections to ensure the derivative is covariant under Lorentz transformations. 
With a metric-geometry described by ( fTSI ) the fermionic terms in the action must be rewritten in terms of the vielbeins and spin 
connections 



24M3' 

= Sl ojs = 0, (16) 



where capital letters from the start of the Latin alphabet have been used to label the vielbeins' Lorentz indices and the Greek 
alphabet characters fi, v are coordinate indicies. Explicitly these become incorporated into the spin-covariant derivative = 
+ (JL>N and gamma matrices T'^ - V'^T'^. With this simplified notation the fermions contribute to the warped space 5-d action 
according to; 

S^^ ^ d^x^{i^T'^DN'V -hx^r''T''Xah'V -l^^'V). (17) 



We use an analysis presented in 02611 to qualitatively demonstrate that if we decompose ^{x, y) according to ([8j, then in warped 
space time there is no mass gap between the zero mode left handed fermion and the continuum of 3-1-1 -dimensional massive 
left and right chiral modes. To make this statement transparent we use separation of variables in the Dirac equation, to obtain 
eigenvalue equations for the extra dimensional profiles. However we now consider the extra dimensional profiles of the massive 
modes, f„u, as well as the zero mode. These satisfy the equations: 

= 6M Vo'i,. - P'hu + 6M'gj<f'jfou « = 

24MSi2eKv)/6M'^„^,. ^ _24M'^/;;^,. + IOMVXl, + [4MV'' -p'^ + 24M''W n>0 ^ ' 

where the generic scalar field coupling term gjcpj = hxii^i'^Xabiu + h^4>(y) has been used to simplify the expression for W - 
{gj(t>jf' - gj<p'j + gj^Pij^ and m„ is the mass of the nth 3-1-1 -dimensional chiral mode. That is: 

iydfti^nLi = m„i//„Ri, ij''d^\j/„Ri = m„ip„u. (19) 
We can solve for the zero mode extra dimensional profile immediately, 

Mi(y)-Noe<'''' (20) 

For the solution presented in (fT4l i there is a normalised 3-1-1 -dimensional left handed zero mode confined to the plane parameter- 
ized by my - In [root (iS/j^M^f^ -i- (6 yf5hx^a.b (r"T*).. {i5aM\ - i6a,b-\)M^ - 5v mjt - IS/z^M^)] provided 6M^h^ - Win > 0. 

For an arbitrary point in the parameter space which will give rise to a different solution to the Einstein-Klein-Gordon equations 
we can still say something about the existence of a confined zero mode fermion. This is because square integrability of (l20l i 
depends on the asymptotic properties of p and (pj. We know that in the case of a scalar field Lagrangian, solving the Einstein- 
Klein-Gordon equations in a vacuum yields a warp factor p (x.\y\, hence we expect p to approach |y| as y — > oo. We also impose 
the boundary conditions Ihriy^-tooX - and //my^-tco^ = ± V, thus any solution to the Einstein-Klein-Gordon equations will have 
the same asymptotic form as ( fT4l i. 

For a continuous, bounded warp factor and Higgs field expressions with the same asymptotic form as ( fT4b there is a delta 
function normalisable zero mode, confined to the y-plane parameterized by the zero of the integrand in ( l20t . provided 6M^h^ - 
Vm > 0. So far our results are consistent with the flat space case. However we must now attend to the massive chiral 3-1-1- 
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dimensional modes. 

Previously we argued that the presence of a mass gap between the zero mode and the lowest discrete mode explained why 
electroweak scale experiments could not detect the tower and continuum modes. It has been shown, by (idn . that in Randall- 
Sundram warped space the continuum modes start from zero mass. To see this we make the substitution fnu ~ ^ p^'' fnu ii^ 
(fTSl l and change variables to conformal coordinates, j: - ePO'^/'^M ^ obtain 



fnLi = mlfnu- (21) 



The change to conformal coordinates, / : y ^ z, is a diffeomorphism from R into a connected segment of the real line so we 
can analyze the potential g ^ ^ function of y and and interpret the results as corresponding to the z-coordinate space 

factored through the mapping : z — > y. This is well defined since is bijective. 

Using our Einstein-Klein-Gordon solutions (fT4l i we can easily see that e"''W/6M ^ asymptotes to zero. Thus we have a 
continuum of delta function normalized modes f„u{z) starting at zero 3-1-1 -dimensional mass, m„, which approach plane waves 
asymptotically. This translates into a continuum of eigenfunctions, f„uiy) Vm„ > 0, which are delta function normalized with 
respect to the weight function g^/'W/^M Davies and George |26] argue that these are precisely the normalization conditions 
required to reduce the kinetic term in the 4-1-1 -dimensional action ( [TtI i to its regular 3-1-1 -dimensional counterpart. We take this as 
the condition for proper normalization of the profile functions. Hence the continuum modes f„Liiy)^/^nLiix) and their counterparts 
fnRi(y)4'nRi{x), which havc analogous conditions omitted here for simplicity, constitute properly normalized solutions to the Dirac 
equation. 

Thus there is a continuum of massive chiral 3+1 -dimensional modes starting from zero mass. It is argued that provided 
< ^ << 1 the potential, g^/'O)/*'*?'^^ will decay slowly to 0. This provides a wide barrier for the low energy, asymptotically 
free, continuum modes to tunnel through. Hence the corresponding wave functions will be heavily suppressed at the position of 
the brane: y — 0. The same behavior is exhibited by the spectrum of Kaluza-Klein modes for the general linearlized fluctuations 
around the metric in the Randall-Sundrum delta function brane case We argue that it is possible for the cross section of any 
process involving interaction between the zero mode and light continuum modes to be imperceptibly low. 

Because the i/'„,l\r,; are not SO(IO) gauge singlets, individual modes will propagate as constituent particles of gauge singlet 
states in the confining bulk. However as argued in section ( llVBb this does not compromise the integrity of our choice of mode 
decomposition for 'V{x,y). Furthermore bound states comprised of low energy, massive, 3 H-1 -dimensional chiral, continuum 
modes will propagate far from the brane, while the zero mode fermion will be trapped around the 3+1 -dimensional topological 
defect. Hence the cross section for processes involving interactions between the zero mode and the low energy continuum modes 
will still be extremely low. 



D. Additional Symmetry Breaking 

The final components of our model are the SU(5)-breaking and electroweak breaking Higgs fields, ((x^^,y) and T](x'',y) respec- 
tively. We introduce these fields now. There are many well documented ways of breaking SU(5)xU(1)a' down to the electroweak 
gauge group and beyond. We choose a specific scenario to illustrate the general process. 
The symmetry breaking pattern: 

SO(IO) D SU(5)xU(l)x 

SU(3)c X SVw(2) X Uf(1) 
— . SU(3)cxU(l)e (22) 

can easily be achieved by using a pair of 16 dimensional representations for SO(IO). The branching rules for this representation 
areiH: 

SO(IO) D SU(5)xU(l)z 

16 10(-l) + 5*(3)+ l(-5) 
SO(IO) D SU(3)cxSU(2)vi.xU(1)fXU(1)x 

16 ^ (3,2)(1)(-1) + (3M)(-4)(-1) + (1,1)(6)(-1) + (3M)(2)(3) + (1,2)(-3)(3) + (1,1)(0)(-5). (23) 

We would like the component of ^(x'^,y) which transforms like (I, l)(0)(-5) under SU(3)cXSU(2)vvxU(l)i'XU(l)A' to condense 
inside the domain wall. The other 15 components of ((x^,y) should not condense. This will ensure that the gauge group is broken 
to the standard model on the brane. 
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Similarly for the field ri{x^,y), we would like the field component which is uncharged under the embedded U(1)q belonging 
to the doublet (1, 2)(-3)(3) to condense inside the domain wall and all other components to not condense. This will implement 
electroweak symmetry breaking on the brane. The most general 4th order SO(IO) invariant potential felt by these Higgs fields is 

= l.ij{Am,h,,h h\ h j + AH2,hi.hj {h\ h >) + .^m [h\ h,) [h\h ^ + Am,h,AiK ^ i^^ 

+AmM.h,h\'^,;^rX^ + AmM.,hh\iT"T^Xahfhj + Am,h„h,h^T^T''Xab(f>hj], (24) 

where the indices /, j run over the set {1,2) with hi = ({x^,y) and = ri(x'',y). 

We factor the A:th entry in the 16-dimensional column vector for each Higgs field into the product of an extra dimensional 
profile function gh„k,n(y) and a complete set of solutions to the 3+1 -dimensional massive Klein-Gordon equations, Ohi,k,n{^^)- 
That is, 

[h[\k = %^ 8h,,k,n(yWh:,k,n(x''), whcrC d^,d^0h„k,n + mX. Sh:,k.n = 0. (25) 

Thus if we let = Uh,,k,hj,p[hj]p + 0(hj), where the indices k and p are being used to label the 16 components of each Higgs 

field, while ; and j still run over the set {1,2) with the same definition as before for hi(x^,y) and h2ix^',y), then we must find the 
solutions to the Euler-Lagrange equations: 

_ d Shnkjy) ^^^^^^^_^-^ ^L—gi^^j^^^{y)0^^j^^^(^x} + Uh_Xhj,pghj,p,n(y)^hj,p,n(x) = m\,.y^^^'''^' gh,,k,niy)ShiXn{x}. (26) 



To get a qualitative idea of the behavior of solutions we consider the simplest case of equation (1261 1. That is we choose the 
coupling constants in equation (l24l i so that Uhi,k,hj,p = when k p or hi + hj. We would like the appropriate symmetry 
breaking components of ri(x'',y) and ^(x'^,y) to have a Kaluza-Klein mode with tachyonic mass. This indicates an instability in 
the background solution h/ = 0. There is sufficient parameter freedom for this to be possible while keeping + Uhi,k,hi,k - 
potential well centered about the coordinate of the domain wall brane. This is necessary to ensure that the 
condensed component of each Higgs field is localised to the domain wall. For the remaining 15 components of each Higgs field 
it is only necessary to ensure that one of these two conditions is voided so that the field does not condense inside the domain 
wall. A thorough analysis would then require us to go back and solve for both the flat space domain wall configuration and the 
Einstein-Klein-Gordon equations consistently with this new background to ensure the backreaction of having additional Higgs 
fields with tachyonic components does not destabilize the brane. This presents a considerable computational task and we do not 
attempt it here. 

The zero mode fermions can now acquire masses through coupling to the SU(5) and electroweak breaking Higgs fields. 
Coupling constants in the effective 3 -t-l -dimensional theory arise from the putative 5-dimensional coupling constant multiplied 
by the overlap integral of the bulk profile functions. This follows from integrating out the y-dependence of the terms in the 5- 
dimensional action. Hence fermions with different bulk profile functions will have different tree level masses. Since the fermion 
profile functions are split along the bulk according to their U{l)x charge the normal SO(IO) tree level mass relations are reduced 
to the less phenomenologically infringing SU(5) mass relations. 



V. SO(IO) ^ SU(3)c X SU(2)z. x SU(2)« x U(1)b-l MODEL MODIFICATIONS 



We would like to engineer a Higgs potential which breaks SO(IO) -> SU(3)c x SU(2)r x SU(2)l x 1](1)b-l directly with the 
background SO(IO) adjoint Higgs. Unfortunately, for our fourth order potential, any Higgs solution respecting this symmetry on 
the domain wall is perturbatively unstable against further condensation. This point is considered in more detail in the Appendix 
(IVIIIb . However it is possible to break SO(IO) down to the left right symmetric model using a sixth order potential: 

+ d(02 _ ^2^2(^2 _ ^2) ^ ^^2^^^2 ^ ^(TrA:^)^' + ^TrAT^ + ^ ArA-^ (27) 

This introduces a plethora of dimension-full coupling constants. However, since we are unable to renormalize 4+1 -dimensional 
Yang-Mills gauge theories, we have assumed there is an inherent cutoff scale beyond which our model is no longer effective, so 
we are at liberty to add sixth-order terms. 
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The Euler-Lagrange equations still exhibit flat space solutions of the form 

(f>{y) = V tanh[my], Xi - Asech[my], 



(28) 



with = -2yU^ + /c +/3V'' and - — ""^ ^K+iw^p"^ ^ "^'^''^ • relabeled the SO(IO) generators so that T\ now refers 

to T\ - -^diag(cr2, 0-2, 0-2, 0, 0). However the conditions are now a lot more convoluted with the solution contingent on: 



(^3 + 2^4 + ^6As)A^ - (48^6 +8/l7V^)A^ + 48/?v'* = 0, 

0, 
0. 



Vln?- + (6/I1 + /I2) - 6/f +6A A^ + A^ - 12y3v'* 



72^ V"* + (6 A + ^7) A"* - 24 fiA^ 



(29) 



A general perturbative linear analysis must now be done numerically, as is explained in the Appendix dVIIII I. However we note 
that for the parameter regime 



2/liA^ + A2A^ - Ik +2^6A^ +^7A^ -4j6v'* > 0, 
(SAsA"* + I4A4A'* + Se^sA"* - 48^6 A^ - l\Ai A^ + 48/3v^) > 0, 



(30) 

it is possible to guarantee that all perturbations about ( 1281 ) are oscillatory and SO(IO) will break stably to SU(3)c x SU(2)£ x 
SU(2)r X U(1)b-l- We cannot predict the behavior of the perturbations outside this parameter regime. A general method for 
identifying the bifurcation points has yet to be determined. The fermions are localised by exactly the same techniques as in the 
SU(5)xU(l)x case. 




FIG. 2: The graph displays the integrands in the exponent of equation JlOt associated with extra dimensional profile functions for the left right 
symmetric model bulk fermions. The graph shown here corresponds to the parameter choices = = ^ and the flat space solution to the 
Euler-Lagrange equations X\ = Asech(ray), if> = Vtanh(mv) with A = V = 1. 3+1 -dimensional left chiral fermions will be confined to the 
hyperplane associated with the zero of their integrand. 



The zero mode fermion profiles are now localised about the x-axis intercepts of Figure|2l 

We can also solve the Einstein-Klein-Gordon equations under these circumstances. We find that our previous ansatz ( fT4l l. 
will satisfy the Einstein-Klein-Gordon equations for the sixth order potential provided we impose the same bulk cosmological 
constant fine tuning and we are in the slice through the parameter regime : 

- 3Af - 12/3 +3A -6k + 12 A^ + 2Ai +2A2 + l2Ai = 0, 

3M^mH V W -I- 6MV^ - 3M^/^ -SM^ySv'* = 0, 

^3 + 6^4 + 36^5 - 12/16 - 2Ai + 3/5 = 0, 

6M^n? + V^n? - 3M^A v'* +3M^A f + 3M^k +6M^/3v'^ = 0, 

9/1 H- 12.16 + 2.I7 - 12yS = 0. (31) 



This symmetry breaking Higgs pattern selects a left right symmetric model on the brane ll28ll . It is subsequently possible to 
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break our left right symmetric model down to the standard model by a technique similar to the one outlined in section ( II V Dl i. 

VI. FURTHER REMARKS 

There are a couple of salient points for 4+1 -dimensional domain wall brane models, based on the Dvali-Shifman mechanism, 
which we have not yet considered. Firstly, since we are constructing an SO(IO) grand unified theory we need the standard model 
coupling constants to unify at some high scale. Secondly, it was mentioned in section dllll l that we need a hierarchical ordering 
of the inherent physical scales for our model to work. Both these points were discussed in the DGV paper lfl2ll . However in 
going from SU(5) to SO(IO) gauge invariance in the bulk we have introduced a new intermediate symmetry breaking scale. For 
this reason we explicitly review both ideas in the context of our model. 

The spectrum of Kaluza-Klein modes will affect standard model gauge coupling constant running. In particular Kaluza- 
Klein modes for 4+1 -dimensional fields carrying different U(l)x charges which belong to the same SO(IO) multiplet in the 
unified SO(IO) gauge theory generally have different masses, i.e. the Kaluza-Klein modes originate in the SO(IO) gauge theory 
as split SO(IO) multiplets. Because the Kaluza-Klein modes form split SO(IO) multiplets, and at specific energy scales only 
certain U(l)z components of the split multiplet will be able to contribute to the beta functions, the comparative running of the 
standard model coupling constants will change. Although coupling constant unification is ruled out for 3+1-dimensional non- 
supersymmetric SO(IO) GUTs, it is still possible for our model. The full calculation would require us to do a phenomenological 
parameter fitting: this will fix the higher mass Kaluza-Klein modes and therefore determine the coupling constant running. We 
have not undertaken this particular task, however we find it tantalizing that the coupling constant may unify at some high scale 
M- 

To maintain the internal consistency of our model, we require a specific ranking of the magnitudes of the inherent physical 
scales. For simplicity we will look specifically at the SO(IO) — > SU(5) model. The concepts do not change for the SO(IO) — > 
SU(3)c X SU(2)l X SU(2)r x U(1)b_l model, hence everything we say here is relevant to this model as well. 

We need to look carefully at 5 specific physical scales: the ultraviolet cut off Auv, the SO(IO) and SU(5) breaking scales on 
the brane Aso(io) ~ [/^{my - 0)]^^^ and Asu(5) ~ (^(i,i)(0)(-5)('M3' = 0))^^-' respectively, the SO(IO) confinement scale in the bulk 
Aconf and the inverse wall width Adw = m. These scales must be ranked according to: 

Adw < Aconf < Asu(5) < Aso(i()) < Auv- (32) 

This hierarchy is necessary for the integrity of the Dvali-Shifman mechanism. Lattice gauge theory simulations |20] tell us 
that the width of the domain wall must be greater than the SO(IO) glueball radius for the Dvali-Shifman mechanism to work. 
This translates into the relation: the inverse wall width is less than the bulk confinement scale. 

The confinement scale in the bulk must be lower than the SO(IO) and SU(5) symmetry breaking scales on the brane. Otherwise 
the non-perturbative bulk SO(IO) confinement dynamics would dominate on the brane and the behavior of the ^ and X fields 
which set these scales would be dictated by their strong SO(IO) interactions. In this case the classical dynamical equations 
we have solved for our background f and X configurations would not be appropriate and our theory would have an internal 
inconsistency. 

Finally our model is an effective field theory which is assumed to be valid only up to the ultraviolet cut off. This establishes 
Auv as the highest energy scale in the theory. 

Amongst the above scales, Aqw, Asu(5), Aso(1{)) and Auv are all functions of the free parameters in the Lagrangian. We argue 
that there is sufficient parameter freedom to rank them according to (|32] |. The bulk SO(IO) confinement scale is putatively to 
be calculated from the value of the dimensionful gauge coupling constant g for the bulk SO(IO) gauge theory. Because this 
gauge theory is not renormalisable in 4+1 -dimensions and we have introduced a UV cut off, Aconf will depend on both Auv and 
g. If we follow the trend from section dllll i and assume that our model has a transition to a confining regime for values of the 
dimensionful gauge coupling constant g, greater than a critical value gc (Auv), then the bulk confinement scale will be set by 
Auv and a g > gc (Auv)- This calculation is beyond the scope of the present work. 

VII. CONCLUSION 

We have identified all the ingredients needed for a 4+1 -dimensional field theoretical model to produce an effective theory of 
3+1-d fields localised to a domain wall brane with SO(IO) gauge invariance in the bulk. We have presented two models, whereby 
the symmetry on the brane is spontaneously broken to give an SU(5) x U(l)x unified model and a left right symmetric model 
respectively. Our models also exhibit localisation of the zero mode linearized fluctuation about the metric. We reiterate that the 
signature points of our model are: 

• All 4 spatial dimensions are infinite. 
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• The brane is supplied by a dynamically generated domain wall. 

• Standard Model fields are localised by their interactions with the domain-wall-engineering Higgs fields. 

• All Standard Model fermions get unified in a single SO(IO) representation in the bulk. However SO(IO) mass relations 
are conspicuously absent from our model. 

Further investigation into the implementation of Dvali-Shifman gauge field localisation arguments in 4-1-1 -dimensions will be 
necessary to establish a rigorous foundation for future models. 

Acknowledgments: We thank R. Davies, D. P. George and N. Pesor for useful discussions and help with computational work. 
JET also thanks H. B. Thompson and Michal Feckan for help with constraining eignevalue spectrums for ODEs. This work was 
supported by the Australian Research Council and the Puzey Bequest to the University of Melbourne. 

VIII. APPENDIX A 

We are interested in the stability of a solution to the Euler-Lagrange equations following from equation 

= aX„b-2iJ^Xab+2ArXab{Xij'^ji) + 2A2{?(')^^^+Kcf>^?(ab. (33) 

We shall call this solution (Xo{y), (po(y))- For the cases examined in this paper we shall only consider Xoiy) for which all the 
nonzero components are in the directions of the elements of the Cartan subalgebra for SO(IO) which we explicitly describe by 
Ci - jdiag(cr2, 0, 0, 0, 0), . . .C5 = 5diag(0, 0, 0, 0, 1x2). To establish the stability of this solution we consider the dynamical 
evolution of a perturbation (d?({x,y, t), S<p{x,y, t)) in a background formed by (^'0(3'), ipoiy))- It is argued that if a perturbation 
is not allowed because of the symmetry of the background domain wall brane configuration or if a perturbation has oscillatory 
time dependence so that ^'^j^fp, 'j^f « 1, Vf > 0, then the domain wall brane formed by (Xo(y),(f>Q(y)) is stable. The 
perturbation 6X(y, t) is allowed to vary independently in the direction of each of the 45 components of the adjoint representation 
for SO(IO). We use the notation 5X{x,y, t) - 6X{x,y, t) + N{x,y, t), where SX{x,y, t) is a matrix containing the perturbations 
along the directions of the Cartan subalgebra generators for SO(IO), and N{x,y, t) contains all the perturbations in the subspace 
of the adjoint representation that is orthogonal to the Cartan subalgebra under the trace operator. Thus N{x,y, t) is a matrix with 
2x2 blocks of zero matrices along the diagonal. 

We will break this argument down by ruling out classes of perturbations. We will start by arguing that perturbations which 
break the 3-dimensional rotational invariance respected by the solution {X^)(y), <poiy)) are not allowed and hence we will suppress 
the coordinate label x and write (Xiy, t), (p(y, t)) = (^'0(3', t) + 6X(y, t) + N(y, t), (poiy, t) + 6(p{y, t)). We then successively rule out 
all perturbations to Xaiy) of the form N(y, t). Finally we will show that all perturbations of the form {6X(y, t), d(p(y, t)) are 
oscillatory in time and hence the background solution (A'()(y), (poiy)) is stable under dynamical evolution. 

Since our solution (Xo(y), <p[)(y)) is invariant under 3-dimensional spatial rotations there is no preferred direction in the hy- 
perplane orthogonal to the y-coordinate. We argue that because the background solution treats all spatial directions orthogonal 
to the bulk coordinate, y, as the same, it does not make sense to say a perturbation has formed along the direction of a specific 
vector in this orthogonal space. This is because all directions are relative to an arbitrary choice of reference coordinate system 
llsoll . Hence we only consider perturbations of the form {SX{y, t) + N(y, t), 5(p{y, t)). 

It is argued in II31I1 that if a perturbation does not change the energy density of our solution {Xoiy), </>o(y)) by introducing a 
term which is linear in the perturbative fields then we can have a stable solution where these perturbations are set to zero in the 
expansion for {X{y, t), (p(y, t)). 

A term which is linear in the perturbative fields would be introduced through an expansion of the SO(IO) Casimir invariants 
appearing in the potential. Hence we consider the form of TrA'^ and TtX'^ for X(y, t) - XQ(y)+6X(y, t) - {XQiy)+6X{y, t)}+N(y, t) 
which we will shorten to A" = [Xq + SX] + N for convenience. For TrA"^ the expression is 

Tr [{Ao + 5X] + N]^ = Tr [{A^ + 2Xo6X + (6X)^} + {Ao + 6X}N + N^] . (34) 

Clearly this Casimir invariant has introduced a term which is linear in the perturbations along the 5-independent Cartan sub- 
algebra direction Ci, . . . , C5. Hence we must examine perturbations along the Cartan subalgebra directions in more detail to 
determine their stability. However for the perturbations along the direction of the other 40 adjoint Higgs field generators, or- 
thogonality of the SO(IO) generators under the matrix bilinear form, trace, means that the linear term disappears. In the above 
equation this is Tr [AqA^] = 0. A similar observation can be made for all the Casimir invariants for SO(IO) because the only pos- 
sible terms in a perturbative expansion which are first order in the perturbations, will be of the form TrA^*^ D Tr Ag*"' (5A -1- A^) , 
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where k eZ^. The exponent, n, for a general Casimir invariants can always be written asn -2k because the antisymmetry of the 
SO(IO) generators implies Tr^V^*^' = . It is easy to check that for any solution Xq written as fields distributed over Ci, . . . , C5, 
<Yg*^"' will also contain nonzero components only in the direction of the SO(IO) Cartan subalgebra generators, Ci, . . . , C5, and 
hence Tr I^^Yq*^ 'A^j - ifsill . This means that the TrA'^ does not contribute any terms to the energy density that are linear in the 
perturbations N(y, t). So a stable solution {Xiy, t), <p(y, f)) can be constructed using only the Cartan subalgebra generators and 
we do not need to worry about linear perturbations along the directions of SO(IO) generators which are orthogonal to the Cartan 
subalgebra Ci , . . . , C5. 

It is left to establish that the perturbations 6X(y, t), along the Cartan subalgebra directions Ci, . . . , C5, have oscilatory time 
dependence under dynamical evolution. For definiteness we choose 6X(y, f) to be the 10x10 matrix: 

6X{y, t) = dXniy, t)Ci + ^^34(3', OCa + ^-^56(3', OCj + 6X7^(y, t)C4 + dX^ioiy, OC5. (35) 

Thus in shorthand our perturbed solution is: 

X(y,t)^Xo(y) + 6X(y,t), cf>(y,t) ^ My) + ^'PiyJ)- (36) 

We substitute this expansion into (|33] | and argue that initially time evolution will be dictated by terms of first order in the small 
quantities S<p{y, f) and 6X{y, t). Discarding all higher order terms we are left with coupled linear homogeneous equations, for 
6(j)(y, t) and the a,bt^ entry of SX{y, t), of the form: 

= D^.^ + K<Po {Xo)ij {5X)j, + {X^^ij (Ao),// + /I {^(Pl - v')) 64,, 

= n6Xat + 4ii {Xo)ab iXQ)ij {6X)ji + {-2^i^ + 2Ai iXo)ij (Xo)ji + ^^,^) 6Xat + 2A2 (6X),i (Xohj {Xo)jb 

+2A2 (Xo)ai (Xo)jb + 2A2 (Xo)ai (Xo)u (SX) j,, + {2k<Po {XqW 54,. (37) 

These equations can be reduced to ordinary differential equations by using Fourier decomposition to factor out the time depen- 
dence of 6(p(y, t) and 6X(y, t). We consider the evolution of a specific normal mode of the coupled system (e"'6X(y), e"'6(p{y)). 

Each mode must independently satisfy the boundary conditions: 5X(y), 6(p{y) — » as y — > ±00 since a global change to 
the profile of either Xiy) or (piy) would require an infinite 3+1 -dimensional energy density and could not be accomplished by 
perturbative effects alone. 

If < for all coupled solutions (e'^'SXiy), e'^'6<p(y)) of the above equations, then the perturbation will be oscillatory and the 
background solution {Xo, <po) will be perturbatively stable over short periods of dynamical evolution. Of course the solutions to 
dJTl i depend on the parameters fp', Ai, A2, A, k. Thus we are looking for relations between the parameters which guarantee < 
for all normal modes. 

Furthermore we are primarily interested in whether the symmetry breaking pattern SO(IO) — > SU(5) x U(l)x is stable. Rather 
than tackling the full problem we will assess when our analytic solution (|4]i is stable, since this will provide an example where 
SO(IO) breaks to SU(5)xU(l)x and stops there. Hence we perturb about {X[)(y), (po(y)) - (Asech(m3'), Vtanh(m3')) and divide our 
problem into two scenarios: either all non-zero components of the superdiagonal of 6X(y) are equal, or two such components, 
lets call them -idXk,k+i(y) ™d -iSX„^„+\(y), are distinct. 

In the first case, there is an unstable perturbation if and only if our solution (HI is not the lowest energy background con- 
figuration for the Higgs fields which breaks SO(IO) to SU(5) x U(l)x and interpolates between the vacua (f> = -\/,X = and 
= +W,X = 0. 

We have numerically solved equation ( |33] ) using a relaxation method for a range of points (fi, Ai, A2, A, k,\/), in the free 
parametric space. We used a trial ansatz X - Xniy)/ VS [Ci + C2 + Ct, + C4 + C5], = ^(y) which automatically selects a 
SU(5) X U(l)x symmetry on the domain wall. The results back up the claim: for a wide range of free parameters the lowest 
energy solution to ( [33T l which breaks SO(IO) — > SU(5) x U(l)x and interpolates between the vacua (p = -\/,X = and 

= +V, A" = has the same qualitative form as (|4]i. 

In the second scenario: there exists k,n e {1, ...9) such that -iSXk^k+iiy) and -i6X„^i,+i{y), are distinct. Therefore we have a 
non-trivial differential equation for dXc = SXk,k+\ - 5A„,„+i. We seek conditions which guarantee aP' <0 for all solutions to this 
equation: 

_ ^ _ ^2m' - j sech2(my) 6X, = {-0? - m^) 6X,. (38) 

We change variables to z = my so that we are working with a dimensionless co-ordinate system. In this co-ordinate system we 
are able to rewrite equation ( |38] ) entirely in terms of a set of parameters with zero mass dimensions {Ai,A2, A, k, V^, A^, w^) — > 
{Aim, A2m, Am, Km, v^/m^,A^/m^, a//m^}. We note that because m has to be real valued to ensure a solitonic background solution 
for (p(y) our condition for perturbative stability, in terms of the new parameter ap-fm^, is aP/nP < . To make it easier to cross 
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reference the results of this stability analysis with section ( IIV Ai we choose not to relabel these new dimensionless parameters. 
Instead of relabeling we point out that stability conditions written in terms of these new parameters will look exactly the same as 
conditions written in terms of the old dimensionful parameters if we set m = 1 . We hope it is clear that this is not a fine tuning 
condition; it is a shift of notation to a convention where all parameters and variables have zero mass dimension. Equation (l38b 
becomes; 



^2-^^yech^z6X,^[-co^-l)6X,. (39) 



Let U - (2 - ^^). Equation (39[ is just a 1-dimensional time independent Schrodinger equation and if we impose the condition 

U > then negative eigenvalues -a/ < will correspond to bound states of the potential - U sech^(z) + 1 . So in what follows we 
are simply determining the conditions necessary for this potential well to be sufficiently shallow so that it does not admit bound 
states. 

We use the following lemma from ll32ll . to conclude that if /I2 > and /laA^ < 10 then there are no nontrivial solutions, to 
equation ( [39] l. with > 0. 



v + (-i + /fsech^f)v' = (40) 



Lemma 1 113211 Let k > 0, A > 0. The equation 



has a bounded solution if and only if there exists an integer M such that 

A = i( V4A- + 1 -4M- If for < M < i( V4/f + 1 - 1) 
or A ^ {{ylAiTTl - AM - if for < M < i( V4^?TT - 3). 

Hence our analytic solutions (|4|i to the Euler-Lagrange equations ( l33T l will maintain a form which breaks SO(IO) — » SU(5) x 
U(l)x under short term dynamical evolution provided /I2 > and /l2A^ < 10. 

If /I2 < then of > Q and there will be a perturbative mode which grows with time. This indicates that our analytical solution 
(|4]l will be unstable under dynamical evolution and it is not possible to stably break SO(IO) to SU(5) x U(l)x with this solution. 
This appears fortuitous as we would expect a cascade effect where the solution for X(y) evolved spontaneously into a form which 
would break SO(IO) to a subgroup of SU(5) x U(l)x on the domain wall. However 6X^ is given by the difference of any two 
arbitrary superdiagonal entries in SXiy). Hence if /I2 < then each nonzero component of 6X will diverge relative to all the 
others. Thus what we will be left with is an SO(IO) adjoint Higgs profile that breaks SO(IO) down to U(l)^ on the domain wall 
brane. Ultimately, mode mixing and other nonlinear effects arising from terms of higher order in SX and S(p will influence the 
long term evolution of the system. Assessing these effect is beyond the scope of our quantitative analysis. We fall back on our 
claim that numerical solutions to equation ( |33T l for X and (p indicate that the analytic solutions (HJi are stable under long term 
dynamical evolution when A2> Q and A2A^ < 10. 

The second half of this paper deals with a sixth order potential for X and (p, dZTl l. where the solution to the Euler-Lagrange 
equations creates a symmetry breaking pattern SO(IO) — » SU(3)c x SU(2)z, x SU(2)r x U(1)b_l. We would like to repeat the 
above techniques with the analytic solution, (l28T l. to these Euler-Lagrange equations. Thereby we will have demonstrated a slice 
of the parameter space where SO(IO) breaks to SU(3)c x SU(2)z, x SU(2)fi x U(1)b_l stably. While most of the results for 
perturbations of the form [dX, 6(f>^ = (dX + N, 6(f>) discussed in the context of the 4th order case can be easily generalized to the 
6th order potential, there will be a slight twist because we can no longer use lemma[T] 

One way of looking at lemma[T]is to make a co-ordinate substitution s = tanhz in equation ( |40] l. This will convert the equation 
( |40] i into a Riemann equation with regular singular points at i = + 1 . The equation is an example of a hypergeometric differential 
equation for which the eigenspectrum is well known I33il . The associated eigenvectors are orthogonal polynomials where the 
inner product is given by 

XI poo 
u(s)l(s)ds^ u(z)l(z) sech^ zdz. (42) 
1 U — cc 

In the case of the 6th order equation if we make the same co-ordinate substitution s - tanhz in equation (l43T l we find that i = ± 1 
are no longer regular singular points and we are no longer able to easily construct a power series solution about these points. 
Hence we do not have polynomial solutions associated with eigenvalues which cause the power series to truncate at a finite 
degree and the previous ideas will not work under these circumstances. 
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We fix this hole in our analysis. Let dX,, - SXjj+i - 6X„^„+i be nontrivial then for m = 1 we have: 

oz^ \ 2 2 I 

1 5h^, ^ IM^, ^ 3^^4 _ ^2 ^2 _ ^2 ^2 ^ sech^z^^,. (43) 
y 48 24 4 2 y 

This is just a 1 -dimensional time independent Schrodinger equation with potential 

U = l+l^A,A^ + y^' - ^2 +AeA^ +^A^ -2ySv4)sech2z + (^A^ + ^A^ + ^A^ - A^ - ^ A^ + ySv^j sech^z. 

(44) 

If we impose the conditions featured in ( [30l l then the potential will be positive definite and the only admissible solutions will be 
bound states with positive energy, hence ap- <Q for all modes, under these conditions. 



[1] T. Kaluza, "On The Problem Of Unity In Physics," Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys. ) 1921 (1921) 966. 

[2] O. Klein, "Quantum theory and five-dimensional theory of relativity," Z. Phys. 37 (1926) 895 [Surveys High Energ. Phys. 5 (1986) 241]. 

[3] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429 (1998) 263 | arXiv:hep-ph/9803315J . 

[4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370 | arXiv:hep-ph/990522]J 

[5] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 4690 l,arXiv:hep-th/9906064| . 

[6] I. Antoniadis, Phys. Lett. B 246, 377 (1990). 

[7] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 257 (1998) |arXiv:hep- ph/9804398 1 . 

[8] M. Visser, Phys. Lett. B 159, 22 (1985) |arXiv:hep-th/99 10093 1. 

[9] V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 125 (1983) 136. 

[10] K. Akama, Lect. Notes Phys. 176, 267 (1982) |arXiv:hep-th/0001 1 13 1. 

[11] G. W. Gibbons and D. L. Wiltshire, Nucl. Phys. B 287, 717 (1987) |arX iv:hep-th/0I09093| . 

[12] R. Davies, D. P George and R. R. Volkas, Phys. Rev. D 77 (2008) 1 2403T|'arXi v: 0705 . 1 5 84 [hep-ph] ] . 

[13] G. R. Dvah and M. A. Shifman, Phys. Lett. B 396 (1997) 64 [Erratum-ibid. B 407 (1997) 452] |arXiv:hep-t h/9612128J . 

[14] E. M. Shin and R. R. Volkas, Phys. Rev D 69, 045010 (2004) |arXiv:hep-ph/0309008 |. ' 

[15] A. Davidson, B. E Toner, R. R. Volkas and K. C. Wall, Phys. Rev D 65, 125013 (2002) [arXiv:hep-th/0202042| . 

[16] D. Curtin, |arkiv :hep-ph/0603217 

[17] G. 't Hooft, Nucl. Phys. B 190, 455 (1981). 

[18] S. Mandelstam, Phys. Rept. 23, 245 (1976). fnlfTil 

[19] N. Arkani-Hamed and M. Schmaltz, "Eield theoretic branes and tachyons of the QCD string," Phys. Lett. B 450 (1999) 92 
|arXiv:hep-th/9812010|. 

[20] M. Laine, H. B. Meyer, K. Rummukainen and M. Shaposhnikov, "Effective gauge theories on domain walls via bulk confinement?," 

IHEP 0404 (2004) 027 |arXiv:hep-ph/0404058 |. 
[21] M. Creutz, "Confinement And The Critical Dimensionality Of Space-Time," Phys. Rev. Lett. 43 (1979) 553 [Erratum-ibid. 43 (1979) 

890]. 

[22] Damien P. George "Domain- Wall Brane Models of an Infinite Extra Dimension" PhD Thesis, University of Melbourne, 2009. 
[23] L. F. Li, 'Group Theory Of The Spontaneously Broken Gauge Symmetries," Phys. Rev. D 9 (1974) 1723. 
[24] V. A. Rubakov Phys. Usp. 44 (2001) 871 [Usp. Eiz. Nauk 171 (2001) 913] |arXiv:hep-ph/0104152J. 

[25] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity, lohn Wiley & Sons, 1972. 
[26] R. Davies and D. P George, Phys. Rev D 76 (2007) 104010 I arX iv: 0705 .1391. [hep-ph] ] . 
[27] R. Slansky, Phys. Rept. 79 (1981) 1. 

[28] G. Senjanovic and R. N. Mohapatra, "Exact Left-Right Symmetry And Spontaneous Violation Of Parity," Phys. Rev D 12 (1975) 1502. 
[29] R. Davies "SU(5) grand unification in an extra dimension" Master of Science Thesis, University of Melbourne, 2007. 
[30] L. Pogosian and T. Vachaspati, Phys. Rev D 62, (2000) 123506 |arXiv:hep-ph/0007045 |. 
[31] L. Pogosian and T. Vachaspati, Phys. Rev D 64 (2001) 105023 |ar Xiv:hep-th/0105128 1. 
[32] K. Yagasaki, Nonlinearity 12 (1999) 799 

[33] A. Nikiforov and V. Uvarov "Special functions of Mathematical Physics" Science, Moscow, 1978. 



